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Abstract 



We explore, within the warped extra dimensional framework, the possibility of finding anti-matter signals 
> ' 

in cosmic rays (CRs) from dark matter (DM) annihilation. We find that exchange of order 100 GeV radion, 
an integral part of this class of models, generically results in a sizable Sommerfeld enhancement of the 



annihilation rate for DM mass at the TeV scale. No ad-hoc dark sector is required to obtain boosted 
annihilation cross sections and hence signals. Such a mild hierarchy between the radion and DM masses 
can be natural due to the pseudo-Goldstone boson nature of the radion. We study the implications of a 
Sommerfeld enhancement specifically in warped grand unified theory (GUT) models, where proton stability 

• rH . 

implies a DM candidate. We show, via partially unified Pati-Salam group, how to incorporate a custodial 
symmetry for Z — > bb into the GUT framework such that a few TeV Kaluza-Klein (KK) mass scale is 
allowed by electroweak precision tests. Among such models, the one with the smallest SO(10) (fully unified) 
representation, with SU (5) hypercharge normalization, allows us to decouple the DM from the electroweak 
gauge bosons. Thus, a correct DM relic density can be obtained and direct detection bounds are satisfied. 
Looking at robust CR observables, we find a possible future signal in the p/p flux ratio consistent with 
current constraints. Using a different choice of representations, we show how to embed in this GUT model 
a similar custodial symmetry for the right handed tau, allowing it to be strongly coupled to KK particles. 
Such a scenario might lead to observed signal in CR positrons; however, the DM candidate in this case can 
not constitute all of the DM in the universe. As an aside and independent of GUT or DM model, the strong 
coupling between KK particles and tau's can lead to striking LHC signals. 
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B The volume factor for antiproton propagation: a diffusion model example 

1 Introduction 

In the last few years a host of experiments have provided us with detailed cosmic ray (CR) data 
in the energy range of 10-1000 GeV [DlllElllElEllTlElEKinilll]. The data is interesting for the 
astrophysics and cosmology communities, enabling them to learn about production and propagation 
of particles in the Galaxy. It is also of great interest for the particle physics community, due to 
the anticipation that annihilation of dark matter (DM), possibly consisting of weakly interacting 
massive particles (WIMPs), would generate an observable signal in the CR data. A lot of model 
building effort has recently been associated with the PAMELA [I] and ATIC/FERMI/HESS [2JEJ 
HI [11] measurements. Probably the main reason for the excitement is due to a rise in the positron 
to the total electron flux ratio (positron fraction) in the 10-100 GeV energy range, as measured 
by PAMELA. The rising positron fraction is in tension with common assumptions regarding the 
production and propagation of CR electrons and positrons in the Galaxy (see e.g. [HI [13], and 
references within). 

The rising positron fraction^, though certainly intriguing, does not necessarily imply an "anomaly" 
with respect to what could be expected from standard astrophysics as follows: The actual positron 
intensity does not exhibit an excess when contrasted with model independent calculations [15] . 
which successfully describe the observed abundance of other secondary CR particles, such as an- 
tiprotons. Moreover, since measurements of unstable CR isotopes can be used to infer the cooling 
suppression of positrons at an energy of around 20 GeV, a theoretical estimate for the corresponding 
positron flux can be derived at that energy [15]. Thus, the combination of the predicted positron 
flux and the available e~ + e + data [2j [31 El [H [9] yields an independent estimate of the back- 
ground positron fraction for this energy range. The authors of [15] have compiled the above data 
and shown that it is, in fact, consistent with the PAMELA measurement, leaving little room for 
an anomaly. It is, therefore, conceivable that the rising positron fraction may just imply that the 
currently fashionable diffusion models for CR propagation in the Galaxy are incorrect. 

Even within simple diffusion models, the PAMELA result has been argued to be compatible 
with secondary positrons, provided that the primary electron spectrum is soft [13] [16]. Along these 
lines there are alternative astrophysical interpretations, wherein the positrons are still of secondary 



See, however, [14] for cautionary notes. 
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origin jlT |, fTS I, [T9] . We note that, at present, all of these astrophysical interpretations require further 
assessment in order to verify the compatibility of the rising positron fraction with the CR nuclei 
and antiproton data. In regards to suggested primary injection mechanisms, pulsars have been put 
forth as astrophysical source candidates (see e.g. |20j ) . and models of DM annihilation or decay 
have been proposed as a particle physics explanation (see e.g. [2"T | |2"2" 1 E ffi I2~4" [ 125])- 

A common feature of the DM annihilation models which can account for observable contributions 
of anti-matter CRs is the presence of a large enhancement ( "boost" ) factor in the annihilation cross 
section. This feature can be traced back to the WMAP data, which fixes the annihilation cross 
section at the cosmological epoch to (av) ~ few 10~ 26 cm 3 s 1 . For DM mass in the TeV range, 
the latter number implies that the positron (or other anti-matter particle) injection rate lies orders 
of magnitude bellow the astrophysical background. One widely studied possibility for obtaining a 
boost factor of the velocity weighted current annihilation cross section relative to its cosmological 
value at freeze-out is the so-called "Sommerfeld Enhancement" (SE) [26], which originates from DM 
particles interacting via a light force carrier. Other forms of enhancement have also been studied 
in the literature (see e.g. [271 1251 [29] ) . 

While not currently necessitated by data, it is still an interesting possibility that the observed 
CR fluxes include a primary component from DM annihilation. Furthermore, in wait of future data 
release by the PAMELA and upcoming missions |30[ 131], it is timely to consider theoretically clean 
observables for indirect detection, such as the antiproton to proton, and positron to antiproton flux 
ratios [T5] . 

In this paper, we explore such robust observables in the future measurements using a well 
motivated theoretical framework, namely, that of a warped extra dimension a la Randall-Sundrum 
model (RSI) [32], but with SM fields propagating in it. One nice feature of the warped extra 
dimension framework in light of indirect astrophysics signal is that there is a natural candidate 
for the force carrier of SE, namely the radion, which is an intrinsic component of the theory^ It 
is the degree of freedom corresponding to fluctuations of the size of the extra dimension in an 
RS-type scenario. Radion mass is in principle a free parameter of the theory, but assuming no 
fine tuning (and Kaluza-Klein (KK) scale of 0(3 TeV)) its mass could vary from 0(100 GeV) - 
in which limit it can be consider as a pseudo Goldstone boson (PGB^fj- all the way to the KK 
scale. The precise radion mass depends on mechanism which stabilizes the distance between the 
branes [34j [35l [36l [37l [38J [39], [40] . Also, the radion coupling to other particles are (roughly) given 



2 Based on AdS/CFT correspondence this nice feature of radion as a mediator of SE is dual to dilaton exchange 
in 4D CFT theories of electroweak symmetry breaking with appropriate DM candidate. Reference [33] considered 
dilaton as messenger between SM and dark sector, but did not study the SE from dilaton exchange. 

3 However, unlike other PGB's, the radion can have sizable non-derivative couplings (required for Sommerfeld 
enhancement) even in the GB limit. 
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by mass of the other particles in units of the KK scale. Hence, for a TeV scale DM, radion coupling 
to DM pair is O (1), and radion mass as large a few hundred GeV can give a significant SE. 

We focus here on a variant of the DM model based on a grand unified theory (GUT) model 
within this framework [4lJ^l where stability of the DM is a spin-off of suppressing proton decay. 
The DM particle in this model is a SM gauge singlet and is a non-standard GUT partner of the 
top quark. We incorporate custodial symmetry protection of Z coupling to bottom quarks [H] 
into the above existing RS-GUT model in order to suppress otherwise large shift in this coupling, 
and construct several models of this type. For simplicity, we mainly focus on partial unification 
based on the Pati-Salam group, which captures the major experimental implication; however, full 
unification is discussed as well. 

We also explore the consequences of implementing a similar custodial symmetry protection of 
Z couplings to right-handed (RH) tau's in order to accommodate the possibility of RH tau's being 
localized near the TeV end of the extra dimension and hence having a large coupling to KK particles 
in this model. In such a scenario, DM annihilation can have a large leptonic branching ratios (BRs) 
via Z' - the extra U(l) of Pati-Salam - exchange^ It is interesting that such large leptonic BRs 
can result in indirect detection in CR positron/electrons. We emphasize that, independent of GUT 
or DM model, such a possibility, in turn, opens up new doors for searching for KK particles (for 
example, KK Z) at the LHC through their decays highly-boosted RH taus, which will be a relatively 
clean signal with negligible Standard Model (SM) background. 

We find, however, that models with significant DM-Z' couplings which allow for such an exciting 
astrophysics phenomenology, in general, yield a too small primordial DM density or are possibly 
in tension with direct detection bounds. Furthermore, this scenario seems to require very large 
representations when fully unified into SO(10) and in any case, it is incompatible with SU(5) 
normalization of hypercharge. Thus, even the SM-level of unification of gauge couplings (which is 
automatic in warped models with SU(5) normalization of hypercharge [15]) is not guaranteed to 
be maintained. 

We hence consider other class of models, which can be fully unified into not-so-large SO (10) rep- 
resentations, and furthermore preserves the SU(5) normalization of hypercharge. Thus, SM level 
of unification of gauge couplings is maintained and even unification with a precision comparable to 
the super symmetric SM one might be possible as in reference |46j . It is quite interesting that this 
model actually predicts vanishing DM-2 7 ' coupling so that the above constraints from relic density 
and direct detection are all satisfied, albeit (as a corollary) not leading to exciting astrophysics sig- 
nals in positron/electron channel. Note, however, that custodial Z — > rf symmetry protection can 

4 Based on above discussion, it is clear that radion mediated SE might also be relevant for other RS-type scenarios 
with DM [421143] . 

5 Recall that the DM is a SM singlet so that KK exchange of SM gauge fields is not allowed at leading order. 
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still be implemented in this GUT model so that that the exciting LHC phenomenology associated 
with tau's is possible. 

The outline of the rest of the paper is as follows. We begin in section [2] with a description 
of the model (with more details in appendix) which is a modified version of the warped extra 
dimensional DM model of references [41] . In section HI we discuss implications for cosmology and 
astrophysics. We explore the SE arising in our framework with a (light) radion. We proceed to 
calculate the DM relic density and direct detection cross-sections. The parameter space compatible 
with WMAP observations and CDMS bounds is delimited. A set of benchmark models within 
this allowed parameter space is defined, in which a large SE factor is a natural consequence of the 
setup. Then, we discuss both particle and astro-physics aspects of DM annihilation. Results of our 
detailed analysis are presented. In section [5] we briefly discuss the radion-related collider signals at 
the LHC with conclusions are drawn in section EJ 

2 The Model 

We first present a review of the general warped extra dimensional framework and then of the DM 
model within it. For a review and further references, see the reference [47] . The reader interested 
only in the particle content of the model and the couplings relevant for signals in cosmic ray 
experiments can skip to tables [U [5] and [6] and the comments listed there. 

2.1 SM fields in bulk of warped extra dimension 

The Randall-Sundrum (RSI) framework consists of a slice of anti-de Sitter space in five dimensions 
(AdSs), where the warped geometry naturally generates the Planck- weak hierarchy as follows [32J. 
The AD graviton, i.e., the zero- mode of the 5D graviton, is automatically localized at one end of the 
extra dimension (hence called the Planck/UV brane). If the Higgs sector is localized at the other 
end (hence called the TeV/IR brane) then the UV cut-off for quantum corrections to the Higgs 
mass can be ~ (TeV), whereas simultaneously the AD gravitational coupling strength being set by 
the usual Planck scale, Mp\ ~ 10 18 GeV. Such a hierarchy of mass scales at the two ends of the extra 
dimension is stable against quantum corrections in the warped geometry, where the effective AD 
mass scale (including UV cut-off) is dependent on position in the extra dimension. Specifically, TeV 
~ JWpi e~ knR , where k is the AdSs curvature scale and R is the proper size of the extra dimension. 
The crucial point is that the required modest size of the radius (in units of the curvature radius), 
i.e., kR ~ 1/-7T log (Mpi/TeV) ~ 10 can be stabilized with only a corresponding modest tuning in 
the fundamental or 5D parameters of the theory [3A\ 140] . Remarkably, the correspondence between 
6 In fact with SM Higgs originating as 5th component of a 5D gauge field (As) it is automatically so [48] . 
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AdS5 and AD conformal field theories (CFT) [39] suggests that the scenario with warped extra 
dimension is dual to the idea of a composite Higgs in AD [501 148] . 

In the original RSI model, it was assumed that the rest of the SM, i.e. gauge and fermion, 
fields are also localized on the TeV brane (just like the Higgs). Such a scenario does not have a 
built-in explanation for the hierarchy between quark and lepton masses and mixing angles (flavor 
hierarchy). In addition, the scenario generically also has a flavor and proton stability problems as 
follows. The (effective) cut-off for the entire SM (i.e., not just the Higgs) is of O (TeV) in this case 
so that the higher-than-dimension-4 SM operators induced by the UV completion of RSI will lead 
to too large flavor changing neutral currents (FCNC's) and too rapid proton decay: recall that 
such operators have to be suppressed by, at least, O (lO 5 ) TeV (if they violate CP in addition) and 
~ 10 15 GeV, respectively, to be consistent with the data. The above argument suggests that similar 
problem would be present for the electroweak (EW) sector, a manifestation of the little hierarchy 
problem. 

2.1.1 Solution to flavor puzzle and problem 

It was realized that with SM fermions propagating in the bulk, i.e., arising as zero-modes of 5D 
fermions, we can account for the flavor hierarchy as well |511 152] . The idea is that the effective 
AD Yukawa couplings of the SM fermions are given by a product of the fundamental 5D Yukawa 
couplings and the overlap of the profiles (of the SM fermions and the Higgs) in the extra dimension. 
Moreover, vastly different profiles in the extra dimension for the SM fermions and hence their 
hierarchical overlaps with Higgs, can be easily obtained by small variations in the 5D fermion mass 
parameters. Thus, hierarchies in the AD Yukawa couplings can be generated without any (large) 
hierarchies in the fundamental 5D parameters (5D Yukawa couplings and 5D mass parameters for 
fermions) . 

As a bonus, the above-mentioned flavor problem is also solved as follows. Based on the above 
discussion, we can see that light SM fermions are chosen to be localized near Planck brane in such 
a way that the effective cut-off for them is 3> TeV. In more detail (this discussion will be useful 
in what follows), the contribution of cut-off effects is actually dominated by near the TeV brane 
(where the effective cut-off is of course of O (TeV)), but the operators are further suppressed by 
profile of the SM fermions near the TeV brane. Since the same profiles dictate the AD Yukawa 
coupling, we see that 4-fermion operators have a coefficient ~ 1/TeV 2 x {AD Yukawa) 2 which is 
sufficient to suppress FCNC's: 

V4 M /TeV 2 _ ipj M /Te\ 2 x profiles at TeV brane 

(SM on TeV brane) ~* ' (SM in bulk) 

~ V4 M /TeV 2 x ( AD Yukawa) 2 (1) 
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As a corollary, the SM gauge fields must also propagate in the bulk (hence the scenario is 
called "SM in the bulk"). Thus, the couplings of SM fermions (with different profiles) to gauge 
KK modes are non-universal, resulting in flavor violation from exchange of these KK modes [53] , 
However, there is a built-in analog of GIM mechanism of the SM in this framework [52] [53] [55] 
which suppresses FCNC's. Namely, the non-universalities in couplings of SM fermions to KK modes 
are of the size of AD Yukawa couplings since KK modes have a similar profile to the SM Higgs, 
i.e., gauge KK modes are localized near the TeV brane. Thus, even though the gauge KK mass 
is of O (TeV), FCNC's from their exchange can be adequately suppressed^ Similarly, the KK 
modes induce effects on electroweak precision tests (EWPT), which can be brought under control 
by suitable imposing custodial symmetries [59| 133]. 

2.1.2 Baryon symmetry 

Satisfying the constraints from non-observation of proton decay requires, however, the mass scale 
of new physics to be generically of O (lO 15 ) GeV so that Yukawa- type suppression of cut-off effects 
on top of O (TeV) scale discussed above is not enough in this case. A simple solution is to impose a 
gaugecjfl baryon-number symmetry, denoted by U(1)b, in the bulk and to break it (arbitrarily) on 
Planck brane so that the "would-be" zero-mode gauge boson is projected out. Thus, proton decay 
operators can originate only on the Planck brane, where they are adequately (i.e., Planck-scale 
which is the cut-off there) suppressed: 



2.2 Dark matter from proton stability in GUT 

Extending the bulk gauge symmetry from the SM to a grand unified theory (GUT) is motivated 
by the resulting SUSY- level precision gauge coupling unification [36] , in addition to an explanation 
of quantized hypercharges of the SM fermions. 

However, the extra gauge bosons in the GUT - for example, X, Y in the case of SU(5) - have 
their KK excitations (with a mass of O (TeV)) localized near the TeV brane (even if their "would- 
be" zero-modes can be decoupled by suitable breaking of the GUT). Hence, if the SM quarks and 
leptons are grand-unified as well, i.e., they arise as zero- modes from the same 5D multiplet in a 
GUT representation, then the, X and Y exotic, gauge KK modes will mediate proton decay with 
only Yukawa suppression [beyond their O (TeV)] mass which is clearly not sufficient. 

7 A residual "little CP problem" [56] is still present [55 1 157] in the above scenario which can be amended by various 
alignment mechanisms [561158] . 

8 Global symmetries are expected to be violated by quantum gravity effects. 



q 3 l/TeV 2 
(SM on TeV brane) 



qH/M 2 Pl . 
(SM in bulk) 



(2) 
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2.2.1 Split multiplets 

The solution is to invoke "split" multiplets, namely, we break the GUT group down to the SM 
by boundary conditions (on the Planck brane so that gauge coupling unification still works). We 
can then choose SM quark and lepton to be zero-modes of two different 5D multiplets in a GUT 
representation. The extra (i.e., "would-be") zero-modes with SM gauge quantum numbers of lepton 
and quark, respectively, from the two 5D multiplets can be projected out by the boundary condition, 
i.e., these fields only have only massive KK excitations. In this way, the X, Y gauge bosons cannot 
couple SM quarks to SM leptons (again such a coupling can only arise if SM quarks and SM leptons 
are contained in the same 5D multiplet): see Fig. [TJ 

However, higher-order effects can "undo" the splitting of quark and lepton multiplets so that 
proton decay can strike again - for example, brane-localized mass terms can mix the (KK) leptons 
from the "quark" multiplet (i.e., which contains a quark zero-mode) with the zero- mode lepton 
from the other (lepton) multiplet and similarly mix (KK) quarks from the lepton multiplet with 
the quark zero- mode from the quark multiplet. In any case, we still have to contend with cut-off 
effects giving proton decay. A simple way out is to impose a U(1)b gauge symmetry in the bulk 
as discussed in the case of non-GUT model. Specifically, the entire 5D quark (lepton) multiplet, 
including the KK leptons (quarks) contained in it, are assigned B = 1/3 (0). 



No mixing: only AS = 1,2... allowed 




No mixing: only AB = 1, 2... allowed 



Figure 1: Split multiplets 

2.2.2 Z3 symmetry 

Unlike in the non-GUT model, a (discrete) subgroup of U(1)b has to be preserved during the 
breaking U(\)b on the Planck brane in order to prevent mixing between the KK leptons from 
quark multiplet with the lepton zero-modes from the lepton multiplet (and thus avoid catastrophic 
proton decay). For example, it is possible to require that the U(1)b symmetry is only broken 
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by scalar fields with integer charges, i.e., only AB = 1,2,... operators are allowed. Thus, the 
above-mentioned mixing of "wrong" (i.e., KK) lepton (or quark) with B = 1/3 (or 0) with correct 
zero- mode with B = (or 1/3) is forbidden (see Fig. [1]), even though 4-fermion proton decay 
operators , albeit safe due to the Planckian suppression, are allowed. 

The crucial observation is that, as a corollary, the GUT partners of the SM quarks and leptons, 
i.e., the (KK) leptons (quarks) from quark (lepton) multiplet cannot decay into purely SM particles 
due to their "exotic" baryon-number assignment. Explicitly, the extra particles in the GUT model 
(including X, Y gauge bosons) are charged under the following Z% symmetry 

$ ^ e 2*i{^-B)Q (3) 

(where a, a are the number of color, anti-color indices on $) whereas the SM particles - having 
correct combination of color and baryon-number are neutral under it. Thus, the lightest Z3 charged 
particle (dubbed "LZP") is stable. 

In references [41] , an 5*0(10) model with canonical representations for SM fermions, i.e., in 16 



was presented. It was shown that SM singlet (RH neutrino) partner of can be the LZP and 
is in fact a WIMP and therefore a good dark matter candidate: a spin-off of suppressing proton 
decay (analogous to imparity in supersymmetry) L_ 

3 Partially & Fully Unified Custodial Models 

In models with the canonical/minimal choice of EW quantum numbers, the shift in Zbb resulting 
from exchange of KK modes is typically (a bit) larger than that allowed by EWPT. This shift 
results in 0(5 TeV) lower bound on the KK scale which implies a rather severe little hierarchy 
problem. A custodial symmetry to protect such a shift in Zbb was proposed in reference [H] which 
requires non-canonical EW quantum numbers. 

Here, we incorporate such a custodial symmetry in the warped GUT DM model of reference |41j . 
presenting several models of this type. For simplicity, we mainly work with partially unified, i.e., 
Pati-Salam, gauge group and comment on full unification into 50(10) on case by case basis. Note 
that there is no proton decay from exchange of X, y-type GUT gauge bosons in Pati-Salam model 
so that there is no motivation for incorporating split multiplets and hence for existence of DM of 
this type in this case. However, we always have full unification into 5O(10), where DM emergence 
is a spin-off of proton stability as mentioned above, at the back of our minds. 

It is interesting that such a symmetry can also be extended to leptons in order to protect the 
shift in Z coupling to leptons. Thus, leptons (in particular, r) can be localized closer to the TeV 



9 The tft multiplet being the one giving the LZP follows from its profile being closest to the TeV brane. 
10 In addition to DM, other GUT partners could also give interesting signals (see e.g. [411 160] ). 
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brane, resulting in larger (than canonical) couplings of gauge KK modes to r. A significant DM 
annihilation to r via exchange of KK gauge bosons, therefore, might be possible, which may be 
relevant for the PAMELA rise (or future signals). As further discussed below this possibility is 
typically in tension with the observed DM relic density and with direct detection limits. 

Finally, although we focus here on models where DM is a SM gauge singlet GUT partner of t R , 
it is worth noting that DM could also be GUT partner of (t, b) l instead, depending on details such 
as the proximity of these profiles to the TeV brane. We will defer study of such a possibility to the 
future. 

3.1 Canonical 

Just to get oriented, the canonical choice for representations under Pati-Salam group, i.e., S , C/(4)c , x 
SU(2) L x SU(2) R are in table [10 Namely, LH SM fermions, i.e., SU{2) L doublet quarks and 
leptons, are SU(2) R singlets with T% R = 0. RH quarks and leptons, i.e., SU(2)l singlets, are 
SU(2) R doublets with T^ R = ±1/2 for RH up quark (or RH neutrino) and RH down quark (or RH 
charged lepton). The SM hypercharge is then given by 



where X are the charges under the non-QCD U(l) generator present in SU(4) C , i.e., SU(4) C ~ 
SU(3) C x U(l) x - We have chosen X = diagv^ (-1/3, -1/3, -1/3,1) when acting on 4 of 
SU(4) C such that the normalization for this generator acting on 4 of SU(A) C is TrX 2 = 1/2. This 
combination of T^ R and X corresponds to the SU(5) normalization of hypercharge when fully unified 
into 50(10). Thus this model (at the least) maintains the SM-level of unification of couplings, even 
in the context of a warped extra dimension. 

Thus, we have the breaking pattern: SU (4) c x SU(2) R — > SU(3) C x U(l)y achieved by boundary 
condition on the Planck brane. The Pati-Salam group is preserved by boundary conditions on 
the TeV brane (of course Higgs VEV breaks SU(2) L x SU(2) R -» SU{2) V ). The gauge field 
that corresponds to the combination of T% R and X which is orthogonal to hypercharge will be 
denoted by Z' . The couplings to Z' are then given by (up to overlap factor denoted below by a) 
(g R / cos 9') (T 3R - Y sin 2 9') , where sin 2 9' = (§c/f ) / (fgf + 9 R ) and g R , 54 are the "4D" couplings 
of SU(4)c and SU(2) R gauge groups, respectively (obviously the normalized U(l)x gauge coupling 
is same as the SU{4) C one). 

Note that, due to Pati-Salam being only partial unification of SM gauge groups, the SU(2) R 
and 5f7(4) c gauge couplings are unrelated so that sin 2 9' is a free parameter. However, it was shown 

Of course, we can invoke split multiplets so that there can be two - one for quarks and one for leptons - multiplets 
of each type in the table. 



Y 




(4) 
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in reference [46J that a 5O(10)-type completion of Pati-Salam, i.e., full unification of SM gauge 
groups, is very well-motivated due to the SUSY- level precision of the gauge coupling unification. 
With this result in mind, we can set 54 = gji to find sin 2 9' = 3/5. 





SU(A) C ~ SU(S) C x U(l) x 


SU(2) L 


S*7(2) L 


LH 


4 ~ 3 1+ li 


2 


1 


RH 


4 ~ 3 1+ li 


1 


2 


H 


1 


2 


2 



Table 1: Canonical representations for SM fermions and Higgs: the subscripts denote the a/8/3 X 
charge. 

3.2 Custodial Pati-Salam model 

As outlined above, we begin by constructing a model with custodial symmetry for Zbb based on 
partial unification, namely, the Pati-Salam gauge group: SU{4)c X SU(2)l x SU(2)r we later 
discuss how to fully unify it. For the implementation of the custodial protection for Zbb coupling, 
the required charges are: 

• T3R = —1/2 for (t,b)L and thus T^r = 0, —1 for tR and 6r to obtain the top and bottom 
massei* 2 ! respectively. 

Thus, we must modify the Pati-Salam representations. Moreover, the SU(2)l and SU(2)r 5D 
gauge couplings must be equal. 

However, the above requirement does not completely fix the model: we first discuss the relevant 
parameters left over below and then describe a variety of models with specific choices of these 
parameters. 

3.2.1 Composite charge leptons 

Once we resort to non-canonical representations, we can choose 

• T3R = for tr (and other RH charged leptons) in order to provide custodial protection for 
its coupling to Z as well. 

In this way, tr can be localized very close to the TeV brani^l i.e., we can contemplate larger 
couplings of KK tr to gauge KK modes (in particular, Z'). Since, via AdS/CFT correspondence, 

12 In the model where top and bottom masses are obtained from the same 5D (t,b)L multiplet. 
13 In order to obtain the charged lepton mass hierarchy, en and (ir might have to be localized farther away from 
the TeV brane than the tr. 
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aeing a composite particle of AD strong dynamics, we will refer to 
Then we must choose T^r = +1/2 for (u,t)l to obtain charged 



such a scenario is dual to tr 
this feature as "composite" 77 
lepton masses. 

One may wonder whether this possibility of having a composite tr is constrained by precision 
tests. For instance, virtual KK Z boson exchange will generate 4-fermion operators involving tr 
and other SM fermions. In our case the dominant constraint comes from (e^f^e) (tj^tr) operator 
since the couplings of KK Z bosons to electrons are vector-like in nature, whereas the tau's are RH 
as discussed above. Using the LEP bounds on such contact interactions from [61J, we find that the 
effective scale suppressing this higher dimension operator should be at least 3 TeV. In our case, the 
KK Z coupling to tr is (roughly) given by ~ gzVknR while the electron coupling is ~ gz/^/kirR 
which gives roughly a coefficient of 1/(4 TeV) 2 for this operator for a 3 TeV KK mass scale, and 
hence is consistent with the bounds. However, with a composite tr, constraints from lepton flavor 
violation might still be an issue which can be addressed by gauging (at the 5D level) the flavor 
symmetries |58j . 



3.2.2 DM couplings to Z' 

Of particular importance are obviously the representation and hence coupling of the DM candidate, 
v 1 . Since the couplings of Z 1 are in general of the same form as canonical model (albeit with a 
different sin#') and DM is a SM gauge singlet (Y = 0), its coupling to Z' is proportional to Tsr, i.e., 
the coupling of v' to Z' is vanishing (non- vanishing) for T^ R = (7^ 0). In the case T^' R = 0, the DM 
coupling to the SM Z (of course induced by higher-order effects) is also custodially protected [34] . 
Obviously, the model's phenomenology differ qualitatively depending on whether DM couples to 
Z' (and Z) or not, so that a crucial choice is 

• + vs. T* R = 0. 

In the following parts we discuss two specific models (two more are given in the appendix) which 
demonstrate the essential differences. Due to the fact that our DM is localized near the TeV brane 
(just like other KK's), a non vanishing T£ R would imply a sizable DM-Z' coupling. This case tends 
to yield a too large an annihilation cross section via Z' exchange into electroweak gauge bosons/top 
quarks and hence typically a too low relic density, unless the DM is of 0(100) GeV in which case 
direct detection from Z exchange becomes a strong constraint. Of course, if in addition tr is 
composite, then the DM annihilation into tr's (which do couple to Z') could be significant which 
could be interesting for indirect cosmic ray positron/electron signal. 



14 Note that the custodial symmetry cannot protect shift in Z coupling to LH charged leptons and LH neutrinos 
simultaneously since we require Tm — T%l for this purpose and LH charged lepton and LH neutrino obviously have 
different T3L, but the same Tsr. 
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As we shall see later, it is quite remarkable that our model with the smallest fully unified, i.e., 
SO(10) representations, actually predict T^ R = 0. Thus, leading to vanishing coupling of DM to 
Z' and SM Z, making it compatible with observed DM density and direct detection bounds. 

Finally, in case where T» R = 0) we also require X u = (7^ 0) in order to obtain Y v = 
(again, in general hypercharge is a combination of T 3R and X, but with a different one than in 
canonical model). 

3.3 Model I (a): T£ R ^ and custodial for leptons 

One possible choice of non-canonical Pati-Salam representations satisfying the above conditions for 
cosmic ray signals in positron/electron is given in table [5J The SM hypercharge is then given by 



Y = T 3R + y/l/6X, (5) 
and the DM and t R arise from a 35 of SU(4). The couplings to Z' are then given by 

9L % (T 3R -Y sin 2 9') 



cos I 



as before, but with sin 2 # 35 = (6gf) / (6gf + g 2 LR ) , instead of the canonical value due to the modified 
combination of T 3R and X entering the hypercharge (note that g R of before is replaced by giR due 
to equality of SU(2)i^i couplings), sin 2 9' is a free parameter at the level of Pati-Salam gauge 
group. We will leave a detailed analysis of completing this model into S'O(10)-type full unification, 
including calculation of the resulting gauge coupling unification in this model, for future work. 
Here, we simply note a few features of a potential unification into SO(10). First, such an extension 
seems to require SO(W) representations larger than 560 [62] . Moreover, even if we find such a 
representation, the normalization of hypercharge above is not the usual SU(5) one so that this 
model does not maintain even the SM-level of unification of couplings. 

However, a loop-level matching of the 5D gauge couplings to the observed QCD and SU{2)l 
ones with the assumption of small tree- level brane kinetic terms gives g^R ~ g± (just like the 
canonical SO(W) case). Based on this observation, we can choose gm ~ 94 (i.e., sin 2 9' ~ 6/7) 
as a "benchmark" value for this Pati-Salam model. It is crucial to realize that the above model 
is just one choice satisfying the conditions of custodial symmetry for the Zbb coupling so that the 
value T3R = — 1 (giving Y = 0) for v' (and similarly the value of sin 2 9' , even with the assumption 
of gLR ~ 54) is not unique: see the model below and the two models in appendix [Al 

3.4 Model II: smallest full unification 

We shall now construct a Pati-Salam model based on the 15 representation of SU(4) and show that 



it is compatible with full unification into 5O(10 

lo Sec also reference 63 . 



:3 



The model also has SU(5) normalization for 
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SU(4) C ~ SU(3) C x U(l) x 


S*7(2) £ 


5£/(2) B 




«JO ~ 08 , I4... 


-1 
1 


Q 
O 


{t,b)L 


35 ~ 38 

3 


2 


2 


TR 


35 ~ 1_4,... 


1 


1 


t) l 


35 ~ 1_4,... 


2 


2 


b R 


35 ~ 38 

3 


1 


3 


H 


1 


2 


2 



Table 2: An example for a model with custodial representations for bi and RH charged leptons, 
with non- vanishing v'v'Z' coupling (see Tab. [5]), the subscripts denote the \/8/3 X charge. 

hypercharge, and it predicts vanishing T^ R and hence DM coupling to Z' jZ . 

The Pati-Salam model is shown in table it can be fully unified into the following SO (10) 
representations: 45 for tR and bR, 120 for (t,b)i and the canonical, i.e., 16 for leptons. So, 
RH charged leptons are not protected by the custodial symmetry, but the model can be modified 
easily to include this feature: for example, LH and RH leptons being (10, 2, 2) and (10, 1, 1) under 
5J7(4) C x SU(2) L x SU(2) R , respectively, which fit into 210 and 120, respectively of 5O(10). 
Moreover, the hypercharge normalization 

Y = T 3R -^X. (6) 

is the same as in SU(5) so that this model maintains SM-level of unification of couplings when 
fully unified into 5O(10). 

3.4.1 Vanishing coupling of Z' to v' pair {T^' R = 0) 

Note that the X-charge of v' vanishes (see Tab. [3]) for this choice of t R representation so that the 
v'v' Z' and v'v'Z couplings vanish. Thus, this case might be uninteresting for indirect searches 
for DM annihilation in cosmic ray positrons/electrons, irrespective of custodial symmetry for RH 
charged leptons - that is why we simply chose the canonical representations for leptons in table [3j 
However, as shown below, it may lead to an observed future signal due to anomalously large 
antiproton flux in the hundreds of GeV region and has the benefit of yielding a correct DM relic 
abundance. And, with the custodial symmetry for RH leptons, LHC signals related to composite 
tr become a possibility. 

The couplings to Z 1 are then given by (qlr/ cos 6') {T 3R — Y sin 2 0') as before, but with sin 2 9' 15 = 
i^dl) I (i#4 + 9lr) as m the- canonical case. 
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SU(4) C ~ SU(3) C x U(l) x 


S*7(2) L 


Sf7(2) B 




lO ~ O -4 , 1q... 


1 


1 




15 ~ 3-4 


2 


2 


TR 


4 ~ l l5 ... 


1 


2 




4 - l lv .. 


2 


1 


b R 


15 ~ 3 4 


1 


3 


H 


1 


2 


2 



Table 3: An example for a model with custodial representations for bi which results in simplest full 
unification. Charged leptons are not protected by the custodial symmetry and the v'v'Z' coupling 
vanishes (see Tab.[6|). The subscripts denote the y^8/3 X charge. 



SM t R , (t,b) L ,TR, fJLR,W, Z,h 



Non-SM 


Comments (quantum numbers) 


v' 


DM: exotic RH v (SM singlet) with 5 = 1/3 


<$> 


radion (scalar with Higgs-like coupling to SM) 


Z' 


extra/non-SM 17(1) in GUT 


X s 


leptoquark GUT gauge boson 



Table 4: Particle content relevant for DM (in-)direct detection. 



3.5 Summary of models characterization 



The relevant particle content and their couplings are summarized in tables H] and [5] ([6]) for the 
partial (fully) unifiable models respectively. A few comments about the tables are in order: 

• u' is the SM singlet (i.e., with quantum numbers of a RH neutrino) GUT partner of t J^l but 
with (exotic) baryon number of 1/3. v' R denotes its RH chirality and has a profile localized 
near the TeV brane (like for any other KK mode), irrespective of bulk mass (c) parameter 
; ^| which dictates the profile of t^. 



for this GUT multiplet 
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Following the notation of references |41j . v' R denotes the Dirac partner {left-handed) of 
Its profile does depend on c for tR in such a way that it moves farther away from the TeV 
brane as tR gets closer to the TeV brane - the v' mass (oc this overlap) decreases in this 
process. 



16 Since, with custodial protection of Zbb coupling, (t, V)l can also be close to the TeV brane, it is possible that the 
LZP comes from this multiplet instead of tR. The analysis for the two cases is similar. 

17 We neglect any GUT breaking here in the 5D fermion mass parameters within a GUT multiplet unlike references 
[41] where small splittings of this type were allowed. 

w v' L was used in references [41! for SU(2)l doublet from {t,b)h multiplet. 
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Coupling 


Value (in units of giR^/kirR) 


Comments 




—a,/ cos -1 9' 


^— v"2 




—a cos -1 9' 




Tri^'Hr 


—iat R cos -1 9' sin 2 9' 


~ 1 


(t,b) Lla Z^(t,b) L 


a mt cos- 1 9' (-i-i sin 2 9') 


a(t,b) L ~ 1 
such that a (ti6)i ~ y ^ ~ 1 


(v,t) l1u Z'»(v,t) l 


a( v>T ) L cos" 1 0' (| + | sin 2 0') 


a (^,r)z, ~ To 


r^a u z' »t r 


a TR cos" 1 6*' sin 2 6>' 


< n 
o TB ~ 1 




a^cos" 1 ^ (-1 + isin 2 9') 


< i 

^6r in 

-K 1U 


Zlong.Z'^h 


2 [Pz long . Ph) 




w, + z' wr 

long. H long. 


az>H 2 [P w + P w - 1 

\ long. long. J 


a^'i? ~ 1 


^bFr^ ( ra dion) 


"V 

"AT ( no gLRVkirR) 


A r = y/6Mpie~ k7rR 



Table 5: Couplings relevant for DM annihilation in model with custodial symmetry for Zbb and 
RH charged leptons, with non- vanishing v'v'Z 1 coupling (see Tab. [2]): value of sin 2 9' is 6/7 and 
note that = 1. 



Coupling 


Value 


Comments 


v'jffuXgtR 


VknR ^a tRV , R 


a t R v' R ~ VOt^ 


v 'r v 'r4> 


'% 


same as in Tab. [5] 



Table 6: Couplings relevant for DM annihilation in simplest fully unifiable custodial case (see 
Tab. [3]): value of sin 2 9' is 3/5, but largely irrelevant for cosmology since T£ R = 0. 

• X s (mostly relevant for the unifiable model with no DM-Z' coupling) and Z' (relevant for 
the partially unified model where DM-Z' coupling controls the resulting relic density) are, 
respectively, the non-abelian and U{\) gauge bosons (beyond gluons) contained in SU (4) c and 
have masses (almost) same as those of KK modes of SM gauge bosons (denoted by Mkk)- 

• Neglecting TeV brane- localized kinetic terms for gauge fields, the couplings can be conve- 
niently expressed (as in the middle column of tables [5] and [6] ) in units of g^oVkirR = gsoVk, 
where g$D is the 5D gauge coupling (of mass dimension —1/2) such that g^r> is the coupling 
of the ("would-be" in some cases) zero- mode (and hence is volume suppressed compared to 
95d)- 

• The custodial symmetry for Z couplings to fermions requires the two SU(2) 5D couplings to 
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be equal, but the SU(4)c coupling is unrelated to it. Hence, there appear two <?4d's in the 
table: giR for the two SU (2) groups and <?4 for 5t7(4) group. 

• <74_d's cannot always be equated to the SM gauge couplings since the relation between the two 
couplings depends on presence of tree-level UV brane kinetic terms and also loop corrections. 
A detailed analysis is left for future work, but here we can choose each of the <?4_d's to 
be (independently) roughly between the SM hypercharge and QCD couplings, i.e., 0.35 ~ 
9lr, 94 ~ 1- 

• The factors a's in middle column of tables [5] and [6] come from overlap of wavefunctions in the 
extra dimension of the involved modes. Specifically, for a coupling of 3 (usual) KK modes 
(which are localized near TeV brane), the overlap gives a ~ 1. Then, each time we replace a 
KK mode by SM/zero-mode we incur a "cost" of ^JasM which is roughly the ratio of profile 
of SM fermion/zero-mode at /near the TeV brane to that of a KK fermion (or equivalently 
the degree of "compositeness" of these SM fermions in the dual CFT description). 

• Similarly, /cry is the degree of compositeness of i.e., the ratio of its profile near the TeV 
brane to that of a usual KK fermion (which is localized near TeV brane). With v' R being 
fully composite (i.e., localized near the TeV brane), the particular appearance of [a~T in the 
table is thus explained. 

• We require yjat R d{t,b) L ~ ^/Ykk such that we can obtain top Yukawa - given by Ykk^<H r a (t,b) L 
- of 1: here, Ykk is the coupling of 2 KK fermions to Higgs and we require it to be smaller 
than ~ 1/7 to allow ~ 3 KK modes in the 5D effective field theory 

• The mixing angle sin 2 9' = (6<?|) / (6g| + g"^f>) appearing in Z' couplings is a free parameter 
(since giR is unrelated to 54), but a "benchmark" value for this mixing angle is 6/7. 

• We use equivalence theorem so that W/Zi ong , is the unphysical Higgs. 

• Finally, the coupling of v' to radion has an additional dependence on c for tn only for the 
case m v > ~ Mj^K/VkirR which occurs for c for ij? ~ —1/2 (in the convention that c = 1/2 is 
flat profile for tp). Since we are most likely not interested in this DM mass region, no factor 
of a is shown here in the coupling of DM to radion. 

4 Implications for Cosmology and Astrophysics 

The potentially light radion, an intrinsic ingredient of the model, has significant implications for 
cosmology and astrophysics. The existence of a light degree of freedom opens the possibility of 
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an enhancement factor in the velocity weighted annihilation cross section, relevant for the current 
epoch, compared to the cosmological value at freeze-out. This effect occurs via the Sommerfeld 
Enhancement (SE) [26]. An enhancement is required in order for annihilation signals to overcome 
astrophysical backgrounds, which would drown those signals for a TeV thermal relic with canonical 
cross section (av) ~ few 10 -26 cm 3 s -1 . 

In section |4"TT1 we explore the SE arising in our framework. Requiring a very large enhancement 
dictates special correlations between model parameters, as well as constrains the radion mass. In 
section 14.21 we proceed to identify the parameter region compatible with direct detection limits and 
with the DM relic density implied by WMAP data. We find that a sizable SE factor is possible, 
and that the model consistent with full unification is viable over a large region of parameter space. 

Indirect detection searches in Galactic cosmic rays, including high energy gamma rays and 
neutrinos as well as antiprotons, provide constraints on the viable magnitude of the SE factor. We 
study those limits in sections 14.3.31 and 14.3.41 For antiproton energies e > 10 GeV, no detailed 
assumptions are required regarding the propagation in the Galaxy. We study possible imprints of 
our model in the high energy antiproton flux, accessible to existing and near future experiments. 
We find that in a sizable fraction of our parameter space (with heavy DM and a PGB radion) a 
p/p future signal is quite generic. 

Regarding CR positrons, as discussed in the introduction, an intriguing hint was reported by the 
PAMELA experiment, suggesting a spectral behavior which can not be easily reconciled with simple 
diffusion models of CR propagation pQ. To our view, this latter observation does not necessitate 
an exotic injection mechanism for the positrons, and we dedicate Sec. 14.3.51 for a discussion of 
this point. Here we comment that our benchmark models which survive the requirements from 
direct detection, provide the correct DM relic density, and adhere with collider and precision test 
constraints, do not exhibit a large enough leptonic vs hadronic branching ratio as required to explain 
the positron fraction rise within the commonly adopted diffusion models. 

4.1 Sommerfeld enhancement with a light radion 

In this section we review the computation of the Sommerfeld enhancement factor, relevant for our 
framework if the radion is much lighter than the dark matter particle [26J. Requiring the maximal 
level of enhancement, SE> 10 4 , implies particular correlations between model parameters. We 
outline these correlations and show, in addition, that lower values of SE~ 10 2 — 10 3 are easily 
accessible. 



The Sommerfeld enhancement due to Yukawa interaction is found by solving the ODE [23 





dx 2 



e 



X 




(7) 
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with 



e„ = - 



£0 



m r 



aM 



a 



A 2 
Air 



and with the boundary conditions 

X(x -> 0) 



, xi x ~^ °°) — * sin(e„x + 5) 



(8) 



0) 



Above, M is the mass of the annihilating particles, v is the velocity of each particle in the center of 
mass (CM) frame, A is the Yukawa coupling and m r is the radion mass. The enhancement factor 
is then given by 



SE 



(x-0) 



(10) 



Using ([8j) we find, for our model, 



a 



6.3 • 10" 



1.2 • 10" 



7.9 • 10~ z — 



v 



MY 
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MY 



(11) 



10 



-2 



To get the effective enhancement one needs to average the SE over the DM velocity distribution, 
which we take as Maxwell-Boltzmann: 

f(v) a v 2 e-'» 2 / 2 * 2 , (12) 

where a is the rms velocity, a = <J J dvf(v)v 2 /3 . Here we use a = 150 kms -1 [641165] . Uncertainties 
of 0(1) associated with the value of the DM velocity distribution could modify some of the details 
of our results, notably when a maximal level of enhancement is considered; yet they would not 
change the overall conclusions nor the detailed results in cases where only moderate enhancement 
levels of order SE ~ 10 2 are discussed. 

The Sommerfeld factor attainable for our model is depicted in Fig. [2j In the left panel, the SE 
is plotted in the (M/A r ,m r /M) plane. Resonance branches cross the (M/A r ,m r /M) plane, with 
enhancement factor SE > 10 4 attainable at the peak of each branch and values of SE ~ 10 3 at peak 
vicinity. The location of the i'th resonance branch in the (M/A r , m r /M) plane follows contours of 
constant values of = e^j , with e^ j ~ 0.6,0.15,0.07, ... arising in the numerical solution to the 
Yukawa problem. Using Eq. (fTTI) we see that the resonance branches correspond to parabolas, 



~M 



a 



MY 
KJ 



(13) 
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where the C^s are constant numbers. Sample values are C\ ~ 0.05, C2 ~ 0.01 for the first (upper) 
two resonance branches. We see that in order to obtain SE > 10 3 , significant correlation is required 
between m r , M and A r . Below we exploit this correlation to extract benchmark model points with 
interesting consequences for indirect signatures in galactic cosmic rays. However, since in our case 
no large SE is required the parameters need not be tuned, as to lie on the resonances in the 
effective potential, hence the sensitivity to small changes in the relevant parameter is probably 
only polynomial. The benchmark models we will consider can easily be located on the right panel 
of Figure [2], in which we plot the SE in the (m r , M) plane for a fixed value of A r = 3 TeV. Direct 
detection constraints (discussed in the next section) are also superimposed on this panel. 
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Figure 2: Left: SE factor projected onto the (m r /M, M/A r ) plane, m r ,M,A r stand for the radion 
mass, the dark matter mass and the scale which suppresses the radion's couplings. Right: SE and 
direct detection bounds, projected onto the (m r , M) plane at fixed A r = 3 TeV. 



Finally, note that the constraint discussed in reference [66j from correlation between Sommerfeld 
enhancement and relic density is not applicable in our case since the relevant particles involved in 
the two processes are different. Furthermore, since our force carrier, the radion, is not ultra-light, 
higher partial waves beyond s-wave are negligible in the current epoch annihilation, and constraints 
due to enhanced DM self-scattering [66j EZ] are easily satisfied. 
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4.2 Dark matter relic density and direct detection limits 
4.2.1 Relic density 

As already anticipated the DM abundance is correlated with the DM-Z' coupling size, in particular 
whether T^ R vanishes or not, we discuss the two cases separately in the following. The analytical 
expressions for the various annihilation cross section can be found in [4T], here we only discuss the 
main qualitative feature of the models parameter space. We compute the annihilation cross-section 
using micrOMEGAs 2.2 [68J for the numerical evaluation of the freeze-out DM abundance (for 
simplicity we have set the KK Z — Z' mixing to zero [41, 69J). 

One important feature of our models is that our DM candidate mass, being the RH top partner, 
is correlated with the localization of (or in 4D dual language, the amount of compositeness of) tn, 
which in turns controls the relic density [41j : 

' 0.65 (c + 1) Mkk if c > -0.25 

todm(c) ~ 



0.83 Jc + l M KK if -0.25 > c> -0.5 (14) 



0.83 y'c 2 - \ M KK exp [k-rrR (c + ±)] if c < -0.5 

where Mkk ~ 2.5k (with k = kexp {—kuR)) is the leading order (++) KK gauge boson mass and 
c stands for the i# bulk mass. For instance for c > —1/4 one finds m u i ~ kiril + c)/2 and for 



-0.4 < c < -1/4, m„/ ~ 2ky/\/2 + C. 

In our calculations, we have neglected (for simplicity) brane localized kinetic term (BKT) for 
bulk fields. BKT's can in principle be used to control the total annihilation cross section and direct 
detection rate as follows [70]. First, BKT's for gauge fields tend to lower the coupling of the lightest 
gauge KK modes to other particles localized near the TeV brane. Such BKT's also lower the gauge 
KK mass relative to the going rate, k, mentioned above. However, electroweak precision tests (in 
particular, the S parameter) put a lower bound of a few TeV on the lightest KK scale which is 
(roughly) independent of the coupling of this KK mode. Combining these two features, we see that 
the annihilation cross-section and similarly direct detection can be reduced by BKT's for gauge 
fields. However, k is then larger than a few TeV which might introduce a severe little hierarchy 
problem into the model. In addition, BKT's for fermions can modify the correlation between DM 
mass and localization of tn and, in turn, some of our conclusions. 

We find that DM annihilation into two radions requires the DM pair (fermion-antifermion) to 
be in a p-wave, and hence is suppressed (see also reference |33jj^l. Thus this channel is not relevant 
for calculation of relic density nor does it get SE. 



19 In general, if the interactions respect parity, then only p-wave annihilation of fermion-antifermion into pair of 
identical scalars is allowed [7T1. 
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Non vanishing DM-2' coupling. For T^ R ^ 0, we find quite generally that there is a tension 
between obtaining the correct relic density and being consistent with bounds from direct detection. 
This is associated with the large v'v' Z' coupling, being enhanced relative to the SM gauge couplings 
by the RS volume factor, \fkirR ~ 6 (or as expected via the AdS/CFT correspondence, being an 
inter-composite coupling j^l. Furthermore, the large T£ R (= 2) for 10 of SU (4) and the fact that 
cos 2 9' is significantly smaller than one (see Tab. [5]), cos 2 6'^ ~ j, for model with 35 of 5^7(4) make 
the rates even larger.. Depending on the mass of u' compared to the intermediate particle, v'v' 
annihilate into the SM particles dominantly through either s-channel (for 2M u i < Mz' ) or t-channel 
annihilation (for 2M v i > Mz 1 )- We find that in the former case the Z' becomes broad enough 
such that resonance enhancement of the cross section strongly suppresses the relic abundance, for 
Mz' ~ 2M v i . For M v i <C Mz> the off-resonance cross section is suppressed by M 2 ; /M| ; . Hence, 
for M u i < Mz'/5 the resulting density is in the right ball park. The case with more massive DM, 
2M u i > Mz> has no kinematical suppression factors and yield a negligible freeze-out density. 

We show in Fig. [3] the resulting J7dm^ 2 for partially unified models [v' GlO as in Tab. [8] (35, 
Tab. [2])] as a function of the DM mass. Curves are shown for Mkk =3,4 TeV, where the green 
curve indicates the corresponding relic density only due to annihilation into the EW sector which is 
rather robust, while the blue curve shows how the density is further suppressed when the coupling 
of Z' to top quark pair is added (we used the canonical choice of sin 8' given by setting g^R = #4, 
which is less robust). The annihilation rate is calculated assuming the mass relation of Eq. (1141) . 
and the smallest possible coupling qlr = 0.35. This choice of gin minimizes the rate, i.e., the 
resulting relic densities can be made much smaller by allowing larger gin, but not much bigger 
(which, as we will show below, induces a strong constraint on these models). Other parameters 
were not varied for concreteness, Mz> = Mkk> A r = Mkk 5 and T^ R = 1 for relic density calculation. 
For the model with 10 of SU(4:), for Mkk = 4 TeV, we see that there is a sizable region of DM 
mass, i.e., below 600 GeV, which gives correct relic density. On the other hand, only the small 
region below 200 GeV works for model with 35 of SU(4). However, as discussed in the following, 
typically both these regions implies a too large rate for the direct detection experiments. 

Vanishing DM-Z' coupling. In this case, since DM coupling to Z' vanishes, the dominant 

annihilation channel is via t-channel X s exchange channel into final state heavy quarks, say tutu 

20 A smaller volume factor would thus result in suppression of the DM annihilation and direct detection rates. For 
example, since the focus here is on unification, one can assume that the UV brane scale is actually the unification 
scale, instead of the canonical choice of Planck scale which gives \fknR ~ 6. However, we have verified that, since 
the unification scale is only two orders of magnitude below the Planck scale, the resulting improvement is only 
incremental. Hence the conclusion about viability of these models is basically unchanged. One can in principle 
consider a much smaller RS volume [421172] . However, then the SM level of unification of gauge couplings and hence 
the motivation for considering a GUT model (and, in turn, the above DM candidate) is lost. 
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(see Tab. [6j). As mentioned, the rate is controlled by the amount of compositeness of the RH tops 
which is also correlated (modulo BKT's for fermions) with the DM mass as in Eq. (I14D . Thus, 
within this case an interesting correlation between the DM mass and the resulting relic abundance 
is obtained. This is also interesting in the context of precision GUT which probably requires 
composite RH top however, the issue of precision custodial unification is beyond the scope 
of this project. We show in Fig. U] the resulting ^dm^ 2 for the simplest fully unifiable model as 
a function of the DM mass. Bands for Mz> = 3, 4 TeV are shown where for each gm is scanned 
over the range gLR = 0.35. .1, both within the favored range, other parameters were not varied for 
concreteness. We see that there is a significantly larger region of the parameter space (than in the 
previous models), i.e., a few 100 GeV to a few TeV, which yields the correct DM abundance. This 
feature is due mainly to absence of Z' exchange in the simplest fully unifiable model. However, the 
more important impact of the absence of DM coupling to Z' is that this model is easily consistent 
with bounds from direct detection experiments (cf. other two models). 



a h 2 vs. M DM for 1 reps, (wilh IV1 KK = 3 and 4 TeV) 




a h 2 vs. M DM for 35 reps, fwifh IV1 KK = 3 and 4 TeV) 
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Figure 3: Relic density, Q-h? vs. the DM mass for the partially unified models with v' G 10 (35) on 
the LHS (RHS). Solid curves correspond to Mkk = 3 TeV and dashed ones to Mkk = 4 TeV and 
9lr = 9i = 0.35. Also shown, as vertical lines, are the constraints from direct detection (purple, 
left most vertical line) and the region (gold) where typically no future p/p can be observed. Direct 
detection bound for Mkk = 3 TeV with 10 of SU (4) is not shown because for this case the entire 
range of DM mass considered here is ruled out by the central value of the direct detection bound, 
while Mkk = 4 TeV with 35 of SU(4) is not shown because, in that case, the central value is 
~ 40 GeV which is below the smallest DM mass shown in the plot, i.e., direct detection is a weak 
constraint in this case. 
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Q h 2 vs. M for 15 reps, (with M = 3 and 4 TeV) 
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Figure 4: Relic density, f2/i 2 vs. the DM mass for the simplest fully unifiable model with v' £ 15. 
The bands correspond to varying the coupling in the favored range gin = 0.35 — 1 and Mkk = 
3, 4 TeV. The points relevant to our two bench mark models ("model L" and "model H") are shown 
as the two light and dark green circles respectively. 



4.2.2 Direct detection limits 

Many experiments are underway currently to directly detect dark matter, and still more are pro- 
posed to improve the sensitivity. In order to ascertain the prospects of directly observing v' in the 
model framework we are considering, we compute the elastic i/'-nucleon cross-section due to the 
t-channel exchange of the radion, which is the most important channel, when radion mass is very 
light. The other important channel is t-channel exchange of the Z, which was computed in [414173] . 

While contributions from radion exchange are generic within our framework the ones induced 
by Z exchange (via Z — Z' or via v' R -KK i/£ mixing) only occur in the partially unified model and 
not in the simplest fully unifiable unifiable model (where the v' coupling is custodially protected). 

The cross-section for Z exchange is (roughly) independent of DM mass, but scales as MZv- 
However, the CDMS bound scales as 1/Mdm (i-e., inverse of number density of DM) hence becomes 
dominant at low masses and tightly constrains the light DM region as shown by the (left-most) 
purple vertical lines of Fig. Several astrophysical unknowns (such as the local DM profile, the 
velocity distribution etc.) are involved in converting the direct detection bound into a constraint on 
a microscopical model parameter space (see e.g. [B]). Nevertheless, for concreteness, taking central 
values seriously, only a very small region survives for Mkk = 4 TeV (none for Mkk = 3 TeV) for 
35 of SU(4), and none for 10 of SU(4). Note that the direct detection bound for Mkk = 3 TeV 
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with 10 of SU(4:) is not shown because for this case the entire range of DM mass considered here 
is ruled out by the central value of the bound. For Mkk = 4TeV with 35 of 5'C/(4) the bound is 
not shown because in that case the central value (DM mass~ 40 GeV) is below the smallest mass 
shown in the plot (i.e., effectively the direct detection bound is weak in this case so that the relic 
density constraint is more important). The reason why the model with DM in 10 of 5L r (4) turns 
out to be more constrained by direct detection than the model with 35 of SU (4) is due to the fact 
that T£ R in the former model is twice as large as in the latter model, while the 6' dependence of the 
effective DM coupling to Z (which controls the interaction with the nuclei) cancels (see table [5]). 
We do not include contribution to direct detection from the radion exchange here since it is highly 
model-dependent, and can be easily made sub-dominant for suitable choice of radion mass and A r . 

However, as mentioned already, DM coupling to Z' vanishes for the simplest unified model of 
15 of SU(4:), which means that t-channel exchange of Z also becomes irrelevant for direct detection 
bounds. Hence, for this model, i-channel exchange of the radion is the single most important 
channel, and direct detection bounds can give information for A r and radion mass. In the CM 
frame, in the non-relativistic limit, the elastic cross-section for radion exchange is approximately 
given as 



where \p' u \ ~ M v iv v i, v v > ~ 10 3 is the DM velocity in the CM frame, mjq ~ 1 GeV is the nucleon 



typical magnitude Atv ~ 10 -6 which includes the radion tree-level coupling to light quarks (u, d, s) 
and gluon, and the heavy-quark-loop two-gluon couplings, with the leading parametric dependence 
A at oc mjv/A r . A sub leading dependence on the mass of radion arises because the radion couplings 
to gauge boson pairs depend on m r . All in all, the model parameters enter the direct detection 
computation in the following way 



For heavy DM, one factor of m 2 N should be replaced by a factor of M 2 , arising from the large 
momentum carried by the heavy u' R and entering the nominator of ()15f) . Finally, note that while 
Eq. (|15p provides a reasonable approximation, useful for obtaining an analytical understanding 
of the parameter dependencies of the direct detection constraints, in practice we incorporate our 
model into the micrOMEGAs [68, 75J package and compute the direct detection bound numerically. 
We find that the numerical results follow the parametric dependence given in Eq. (|16p rather well. 




(15) 



mass, \n/\/2 is the effective rNN coupling, and t is the Mandelstam variable, which can be 
ignored compared with m 2 in the radion propagator. For the radion- nucleon coupling we find the 




(16) 
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Our results are illustrated in Fig. [21 where the direct detection constraints are superimposed 
on top of the SE factor. Taking the direct detection limits on face value, we find that a very light 
radion of m r < 20 GeV is already excluded by both the CDMS and Xenon experiments. The 
CDMS limit disfavor m r of up to ~ 40 GeV for DM mass as high as 3 TeV. The entire region of the 
remaining parameter space, where our analysis is valid, will be probed by upcoming experiments, 
such as Super-CDMS and Xenon 1-ton |76j . 

4.3 Indirect detection, simplest fully unifiable model 

In the following sections we evaluate the implications of our framework to various CR species. 
As we have discussed above, models where the DM is not custodially protected are in tension 
with direct detection experiments or lead to too low relic density. We therefore focus on the 
unifiable model where the DM— Z' coupling vanishes. To facilitate the discussion we introduce two 
benchmark models and study the resulting CR injection spectra and rates. We then move on to 
signatures in photons and neutrinos. High energy photon and neutrino observations constrain the 
DM annihilation cross section, weighted by the integral of the DM number density-squared along 
the line of sight of the experiment. 

Proceeding to antiprotons, we note that the astrophysical background is constrained by existing 
CR data. Subject to few general assumptions, the effect of propagation in the Galaxy can be 
accounted for at the cost of introducing a single additional fuzz factor. The antiproton analysis 
is, in this sense, as predictive as the analysis of photon signals where the analogous fuzz factor is 
contained in the line of sight integral. Lastly, we turn to the more involved case of positron signals 
and briefly discuss the injection rate of e + /p for which the astrophysical background is somewhat 
easier to interpret. 

4.3.1 Benchmark models and CR injection spectrum 

Following the discussion of the relic density, Sommerfeld enhancement and direct detection bounds, 
we focus here on two viable benchmark model points characterized by different values of DM and 
radion masses which result in turn with different annihilation spectra. We keep fixed the value of 
the Z' mass, mz> = 3 TeV. The benchmark models are defined as follows 

Model L: M = 600 GeV, m r > 40 GeV, which corresponds to the LH circle on Fig. UJ In 
principle one can obtained a sizable SE while decreasing A#, however, in this case we find 
tension with direct detection bounds (from radion exchange). 

Model H: M = 2400 GeV, m r = 0(100) GeV, which corresponds to the RH circle on Fig. H In 
this case theres a wide range of radion masses and corresponding Aj? which yield a sizable SE 
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and consistent with direct experiment. 



In both cases the annihilation is dominated by v'v' — * tutu via t-channel X s exchange, the couplings 
are given in table [6] and according to Eq. (|14p which link the top compositeness with the DM mass. 

The CR injection spectra of stable final states are plotted in Figure for the various benchmark 
points. These spectra, together with the DM mass and Sommerfeld enhancement factor serve as 
the particle physics input required for the calculation of indirect detection signals. 




1 O 1 1 9 1 

10 10 10 10 10 10 

E [GeV] E [GeV] 



Figure 5: Decayed final state annihilation spectra for the two benchmark models. 



4.3.2 CR production rate 

The production rate of a cosmic ray specie a due to annihilation of Dirac fermion DM at a given 
spatial position r in the Galaxy is given by 

< 17) 
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Here n{f) is the total DM number density (particle+antiparticle) and d(7V ( DM J^ M ^ a ) \ s differen- 
tial velocity-weighted annihilation cross section for the production of the specie a. It is convenient 
to work with dimensionless quantities, 

E M av , n(r) 

GeV TeV 6-10 Zb cm. 6 s 1 n(r so i) 

where M is the DM mass, av is the total velocity-weighted annihilation cross section, f so \ ~ 8.5 kpc 
is the distance between the solar system and the Galactic center and n(f so \) = 0.3cm _3 GeV/M 
is the DM number density in the local halo. For the local halo mass density, we adopt a value of 
Adm(Jsoi) = 0.3 GeV cm -3 . Order one deviations for this number are possible both on average and 
due to local clumps, and go through to the computed CR flux. With the definitions (fTBj) . the CR 
production rate can be written as 

<2 Q ,D M ( e > f ) = Q a , D M^^ so \) x n 2 (r) , (19) 

with the local injection rate 

Q Q ,DM(e,r so l) = 1.3 • 10~ 33 || ^ cm" 3 s" 1 GeV" 1 , (20) 

and where is the differential number of stable final state particles of specie a emitted per 
annihilation event. In writing Eq. (|19p we have neglected the spatial dependence in the Sommerfeld 
enhancement [77] . As, in this paper, we do not attempt to provide a detailed description of the 
spatial features of the DM annihilation signal, we neglect this possible complication throughout the 
discussion. 

The rate of DM annihilation is proportional to the number density squared, and so the results, 
in particular as concerns photon and neutrino flux from the Galactic center region, depend on the 
assumed profile. The latest N-body simulations, including only DM and no baryons, point to DM 
halo profiles with a cusped central region. However, the inner zone of a few hundreds of parsecs 
from the center remains uncertain. In addition, the effect of baryons may be significant at the 
central region and its impact on the DM distribution is far from understood. Baryons were argued 
to either increase the inner cusp, or actually smooth it out, resulting with a cored profile \78\ 179]. 
In this work we analyze both cusped and cored DM halo profiles. The examples we consider are 
the cusped NFW [80] and the cored isothermal sphere [8l] (denoted below by ISO). We do not 
attach special significance to any particular profile but rather lay out the consequences of each 
case regarding indirect detection prospects for our framework. The radial dependence of the halo 
distributions is, 

p(r) _r so \ ( 1 + r sol /r s x 1 



NFW: =— 7 , r s = 20kpc, 

p{r so i) r V l + r/r s J ^ 

iso: J^L = i±r|M, r , = 5kpc . 

p[r so \) 1 + r z /rj 
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4.3.3 Photons and neutrinos 



The incoming flux of photons or neutrinos per unit solid angle is obtained by integrating the 
production rate along the line of sight in a given direction Q in the sky, 



Gamma ray observatories report limits on e), which is defined by averaging (|22p over the 
acceptance Afi of the experiment, 



The observed photon flux depends on the local injection rate, up to a single over all model dependent 
factor given by the line of sight integral, which encodes the DM distribution. 

We derive gamma ray-based model constraints from the following data sets, provided by the 
HESS imaging air Cherenkov detector: 

• HESS observations of the Galactic Center (GC) [52] : 

The GC data set corresponds to the inner 0.1° of the GC gamma-ray source, HESS J1745- 
290. The energy range was £L > 160 GeV. Considering the uncertainties involved in the 
calculation, we find it sufficient for our purpose to use the power law fit reported by the 
HESS collaboration, j oc £T r with T = 2.25 ± 0.04(stat) ± O.lO(syst). The normaliza- 
tion is defined from the reported value of the integrated flux above 1 TeV, J dEjAQ = 



[1.87 ± O.lO(stat) ± 0.30(syst)] • 10~ 12 cm -2 s" 1 . For definiteness, we use the central values 
for the reported flux and impose that the photon flux resulting from DM annihilation does 
not exceed it in order to derive constraints on the model. 

• HESS observations of the Galactic Ridge area (GR) [83] : 

The GR data set corresponds to an observation of the rectangular angular patch \l\ < 
0.8°, |6| < 0.3°, from which the spectral components of the sources HESS J1745-290 and 
GO. 9+0.1 were subtracted. The energy range was > 170 GeV. We use the power law fit 
reported by the HESS collaboration, ] = k (^)~ r with T = 2.29 ± 0.07(stat) ± 0.20(syst) 
and the normalization k = [1.73 ± 0.13(stat) ± 0.35(syst)] • 10~ 8 TeV -1 cm~ 2 s _1 sr -1 . We 
use the central values for the reported flux and impose that the photon flux resulting from 
DM annihilation does not exceed it in order to derive constraints on the model. 

see [84] for the details of the source GO. 9+0.1 




(r). 



(22) 




(23) 



ITcV 



30 



Data from the Fermi-LAT satellite-borne detector has recently become available. We analyzed the 
preliminary results presented in [85] for the Galactic center region. This data constrains the lower 
energy part of the spectrum and, for model L with a cuspy DM profile, is competitive with the 
HESS data. 

The situation is illustrated in Fig. [UJ in which we plot the GC data set of FERMI and HESS vs 
model signals, evaluated with an NFW DM halo profile and the maximal Sommerfeld factor allowed 
by the GC data set. (note that, for model H, the HESS GR data set is in fact more constraining 
and a value of SE = 180, used in the figure for illustration, is excluded.) 

Limits on the neutrino flux arise from measurements of the neutrino-induced muon flux in 
neutrino detectors. For DM mass M, the flux of muons at the detector is given by 



/ defj, / de v n N j Vfl (£l, e)AO 



da vN da pN 



+ 



Here j Uii (£l,e) is the muon-neutrino flux at the earth, which equals the anti-neutrino flux in the 
case of DM annihilation and is obtained from Eq. (|23|) (We use Tri-Bimaximal neutrino mixing for 
definiteness) . The differential cross sections are given by 

da uN _ 2G 2 F m N 



de^ 7r 



2 



(25) 



with a\ ~ 0.2, a% ~ 0.05 for neutrino-nucleon CC scattering and the same with a\ t 2 interchanged 
for the antineutrino-nucleon case. The muon range in the rock beneath the detector is 

L(e„,e th ) = ^-ln( a ^ + ^ e A , ( 26 ) 



where eth = 1.6 GeV is the threshold energy for detection, and a M « 2 • 10~ 3 GeVcm 2 g _1 and 
(3^ ~ 3 • 10 -6 cm 2 g _1 are the muon energy loss coefficients. For the target material we consider a 
nucleon mass rh-N = m p, and the nucleon number density is given by = p/ffiN- 

We derive neutrino-based model constraints from the upper limits on the upward through-going 
muon flux, measured at Super Kamiokande (SK) [86]. We use the 95%CL limits quoted in [87], 
where the line of sight integrals were also given for angular acceptances of 3° — 30° and various DM 
halo profiles. 

In Table [7] we summarize the photon and neutrino constraints. Regarding observations of the 
Galactic center region, the line of sight integral depends on the assumed DM halo profile as well 
as the angular resolution of the experiment. Small changes in the halo profile around the poorly- 
known central regions of the Galaxy result with significant variations in the predicted flux [88j . For 
the cored profile, large cancellations can occur due to background subtraction, and the resulting 
bound becomes weak [87| in comparison with antiproton and neutrino constraints. In this case, for 
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Table 7: Upper bounds on the Sommerfeld enhancement factor SE, resulting from HESS and 
FERMI 7 and SK v constraints. Square brackets refer to optimal background subtraction with 
the HESS resolution. For v we report the result corresponding with the most constraining opening 
angle for the SK analysis. In case the bounds are weaker than 10 4 , we keep only the order of 
magnitude. We also quote the analyses of Fermi constraints, provided in |89l ,90]; see text for 
details. 
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the HESS analysis we report the bound without accounting for these cancellations (given inside 
square brackets in Table [7j), such that the optimal performance can be assessed. 

Coming back to the Fermi data, we note that a part of the power of these measurements lies in 
the complete coverage of the sky. As a result, strong constraints can be derived also for cored DM 
profiles, which were effectively unconstrained by earlier measurements. A more complete treatment 
of the new Fermi data, which included the same final state annihilation products as in our model, 
was very recently provided in |89[ [90] . The analysis of these references is in good agreement with 
ours for observations of the Galactic center, but as expected, it presents much stronger bounds for 
the cored profile. In particular, according to |89[ 190]. the SE for our model L(H) cannot exceed 
~ 500(H) 3 ) in the ISO profile scenario. For an NFW profile, the SE for model L(H) is limited below 
~ 50(100). 

Putting all together and including the results of [89, [90], we find that the neutrino bounds are 
sub-dominant in comparison with the new photon data, for any DM profile. Finally, note that with 
a realistic treatment of the backgrounds, the bounds we apply are likely to tighten by a factor of 
at least a few, implying that the SE factor for our models will probably be limited to a few tens 
(hundreds) in case a cusped (cored) DM profile is adopted. As we show below, such value of the 
SE is still sufficient to produce interesting antiproton signatures. 

4.3.4 Antiprotons 

The PAMELA experiment has recently measured the antiproton to proton flux ratio [6]. The 
reported antiproton fraction does not show deviations from the expected result, based on secondary 
production by pp and spallation interactions of primary CRs with inter-stellar medium (ISM). 
Nevertheless, DM annihilation can contribute a primary component to the CR antiproton flux, 
with a production rate density given by Eq. (|19p . This contribution must be small at the currently 
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explored energies, but could in principle reveal a peak at > 100 GeV energies, soon to be measured 
by the PAMELA and (hopefully) AMS02 [30] experiments. 

Cosmic ray antiprotons have long been considered as a good channel to detect exotic sources 
(see e.g. |91j). At high energies (E > 10 GeV) the background can be determined from the CR nuclei 
data [92] [Td] , leaving significant predictive power. As concerns the DM contribution, analyses 
in the literature were based on detailed propagation models. Such propagation models typically 
include additional free parameters which reduce the predictive power of the analysis. Here we show 
that the antiproton flux can be computed in a model independent manner, at the cost of introducing 
one free parameter to the calculation. This parameter is an energy independent effective volume 
factor, encoding the different spatial extensions of the DM and the spallation sources. The fact that 
only one free parameter is introduced makes the antiproton channel as predictive as the photon 
and neutrino channels, for which propagation in the Galaxy is trivial. 

The approach we adopt is based on the fact that high energy antiprotons above a few GeV 
suffer only small energy losses as they travel through the Galaxy, and on the fact that the secondary 
antiproton flux up to E ~ 300 GeV can be computed in a model independent manner, based on the 
existing CR nuclei data [92j [931 E] • To proceed, we need the following ingredients. 

The first ingredient concerns the propagation in the Galaxy. Define the quantity G(e, es; r, fs) > 
encoding the propagation of CR antiprotons in the Galaxy as follows 

np(c,r) = J d 3r s J de s Qp(es,rs)G(e,es;r,r s ) , (27) 

where Qp(es,fs) is the injection rate density at energy eg at the point fs , the spatial integral 
contains the confinement volume of Galactic CRs and np(e, f) is the antiproton density at some 
point f in the Galaxy ^|. The negligible energy change of antiprotons above 10 GeV implies that 
G(e,es;f,fs) oc <5(e — es) . We now make the assumption that G is separable, i.e. that 

G(e, e s ; f, f s ) = 6(e - e s )g(e)G(f, f s ) ■ (28) 

The second ingredient concerns the secondary source spectrum. We assume that the injection 
spectrum (not rate) of secondary antiprotons has a homogeneous distribution in the Galaxy. In 
practice this assumption amounts to demanding that spatial variations in the spectrum of primary 
CRs are small, at least in the regions from which most of the secondary antiprotons observed locally 
are generated. Under this assumption the secondary source term is separable, 

Qp,sec(£)^) — Qp,sec '"sol) X Qsec(T') , (^9) 



22 For simplicity we did not introduce a time label in Eq. (I27|l . which could account for deviations from steady state. 
Provided that the explicit assumptions we make hold, adding time dependence to the problem would not change our 
basic result. 
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where Qp, se c(e, ^*soi) is the local secondary injection rate. 

Under the above assumptions, the local density ratio between the primary and secondary com- 
ponents takes the form: 

np,DM(e,r so i) Q p , DM (e,r so i) 
— — , ^ , = Jv-fi , j r , 

"pjSeci^) ' solj *Vp,sec l^C) '"sol J 

with the energy independent volume factor, 

/ d 3 rq DM (r)G(r soh r) 
fV = J d 3 rq sec {f)G(r soh r) ' (31) 

For a DM annihilation source, we have qDMif) = w„(r) . We can write the antiproton to proton 
flux ratio as follows, 



Jp(e,r so i) _ Jp{e,r so \) _ 

X 



l + /y 



Qp.DM ( £ > r sol) 
Qp.BCC ( £ 5 r S0l) 



(32) 



«/ P (e,r*soi) V^p(e,r so i)/ sec 

The first factor on the right hand side is the secondary antiproton to the primary proton flux 
ratio. This quantity is constrained by the B/C data, leaving no free parameters. We conclude that 
under some general assumptions, the antiproton to proton flux ratio including a DM contribution 
can be computed based on the relatively well constrained local injection rates and only one addi- 
tional parameter, fy, encapsulating all the details of the propagation. A naive estimate suggests 
fv ~ L/h ~ 10 — 100 , where L ~ 1 — 10 kpc is the assumed half width of the CR propagation 
volume and h ~ 100 pc is the half width of the Galaxy gaseous disc. 

The class of models for which Eq. (1321) holds includes the disc+halo diffusion model with a 
homogeneous diffusion coefficient [94, 95j^j. In Appendix[B]we use this model as a concrete example, 
deriving the precise realization of Eq. (132D . We find, as expected, fv in the range ~ 10 — 100, 
depending mainly on the size of the CR confinement halo with an order one correction depending 
on the DM distribution. 

In Fig. [7]we plot the antiproton to proton flux ratio with a DM component, corresponding to our 
benchmark models. The curves including DM contribution are obtained by suppressing the pure 
background term to 75% of its central value (we find that a similar suppression also best describes 
the data with only the background component), and boosting the DM term by the factor SE xfy, 
indicated in the plot. The shaded region denotes an 40% uncertainty estimate for the background 
calculation [93]. Data points are taken from published [6] and preliminary [97J PAMELA data. 
As illustrated in Fig. [7J a future p signal can arise for ttidm ^ ITeV, with SE x/y > 10 3 . As a 
volume factor fv > 10 is envisioned, the requirement on the Sommerfeld factor is SE > 100, easy 
to obtain in our model with a 100 GeV radion. The TeV scale for the DM mass, roughly above 



23 Of course, the class of models for which Eq. (|32p holds, includes also the well known leaky box model [96]. 
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which the resulting antiproton feature can be pushed higher than existing constraints to provide a 
future signal, is indicated by the vertical golden line of Fig. HI 

Concerning the astrophysical background calculation depicted in Fig. [7J a comment is in order. 
Extending the background prediction all the way to E ~ TeV requires extrapolation of the CR 
grammage (provided in Appendix [Bl Eq. (f3"T|) ) beyond the 200-300 GeV, where reliable data exists 
[98]. While there are indications that the grammage used in Fig. [7J persists to TeV energies |31j . 
the issue is not currently settled [99]. We anticipate that with improved compositional CR data, 
extending to TeV energies, an updated model independent prediction for the p/p ratio will become 
directly available along the same lines described above [92], [93J [15] . Eq. (i32|) will then become useful 
up to TeV energies. 

4.3.5 Electrons, positrons and the positron to antiproton flux ratio 

Recently the PAMELA collaboration have reported a rise in the positron to electron plus positron 
fraction [fj, beginning at E~ 10 GeV. The reported rise have induced numerous publications, 
suggesting an explanation in terms of DM annihilations or decay. However, before examining exotic 
contributions it is necessary to first understand the astrophysical background, which is harder to 
constrain than in the antiproton example. 

In fact, since secondary positrons are produced by pp and spallation interactions, just like 
antiprotons, an upper bound to the positron flux can be obtained model independently, based 
on the measured CR grammage [15] . Contrasted with the data, this calculation reveals that the 
rising positron fraction is not accompanied by any actual positron excess with respect to the model 
independent upper bound. One is forced to the conclusion that the rising positron fraction most 
likely corresponds to an unexpected spectral behavior of the suppression due to propagation energy 
losses, denoted here by f s . Using the total e + + e~ measurements [7] [5] [2] [3] in conjunction with 
the PAMELA data, one finds f s « 0.3 at E 10 GeV, rising to / s « 1 at the highest data 
bin E ~ 80 GeV. At E < 40 GeV, the suppression of the positron fraction can also be compared 
with the measured suppression due to decay of the flux of radioactive unstable CR isotopes, such 
as 10 Be, 26 Al and 36 C1 |100[ 195]. In particular, measurements of the (purely secondary) decaying 
charge to decayed charge Be/B extend to rigidity of ~ 40 GV [100] . These measurements suggest a 
value of f s ~ 0.3 for positron energy ~ 20 GeV, in agreement with the actual result and in support 
of the secondary origin of the detected positrons. 

To summarize, it is our view that the rising positron fraction does not consist an evidence for 
exotic components in the positron flux, simply because there does not seem to be any positron excess 
- merely an intriguing suppression pattern. If, however, future data release by the PAMELA mission 
or other experiments [30] will establish that the rising behavior persists to e + /(e + + e — ) > 0.2 
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Figure 7: The antiproton to proton flux ratio with a DM component, corresponding to our bench- 
mark models. The curves including a DM contribution are obtained by suppressing the background 
prediction to 75% of its central value, and boosting the local DM injection rate by the factor SE x fy. 
This factor encodes the combination of propagation, via the volume factor fy, and of the Sommer- 
feld enhancement SE. The shaded region indicates an 40% uncertainty estimate for the background 
calculation. Data points are from published and preliminary PAMELA data. 
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beyond ~ 100 GeV, a positron excess will indeed be implied, necessitating a primary source. 

Both of our benchmark models do not produce a hard lepton flux, hence, no anomaly is predicted 
in leptonic channels. This conclusion is supported by the expectation that, due to radiative losses, 
positrons do not experience the volume enhancement relevant for antiprotons. However, as argued 
above, since the background distribution (both of primary electrons and secondary electrons and 
positrons) is largely unknown, we proceed to discuss and analyze, in the following part, a more 
robust observable related to secondary to secondary flux ratio. 

In terms of theoretical uncertainties, the positron to antiproton flux ratio is a clean discriminator 
between a secondary astrophysical production mechanism to any other hypothetical source. The 
reason for this is that secondary antiprotons and positrons are produced by the same mechanism, 
namely pp and spallation interactions of primary CRs with ISM. The relative amount of positrons 
and of antiprotons injected at a given energy depends on the corresponding branching ratios and, to 
a lesser extent, on the spectrum and composition of the primary CRs and the ISM. An examination 
of the dependence of the positron to antiproton ratio on the spectrum of primaries was carried out 
in [15], where this dependence was found to be very mild. Since at high energies energy losses 
affect only the positrons and act to suppress the observed flux, and since in the absence of losses 
high energy positrons and antiprotons would propagate in a similar way, the positron to antiproton 
injection rate ratio forms a robust upper bound on the corresponding flux ratio, relatively immune 
to propagation details. 

The spectrum of final state products in DM annihilation may deviate significantly from the 
corresponding branching ratios in pp collisions. The existence of a DM component in the CR 
antimatter flux can therefore be searched for in the positron to antiproton flux ratio. Finding this 
ratio above the standard prediction (based essentially on the branching ratios in pp collisions, with 
mild compositional corrections) will provide strong motivation for an exotic contribution. 

In figure [8] we plot the positron to antiproton production rate ratio for our benchmark models, 
including the background, compared with the prediction for pp collisions which can be regarded 
as an upper bound for the background result. In both models, the ratio lies very close to the 
astrophysical background. The conclusion is that it is unlikely, yet not inconceivable, that our 
models would lead to an excess in leptonic CR signals. 

5 Radion Collider Phenomenology 

For the region of our model parameter space, where a possible CR signal in DM annihilation is 
obtained, a lightish radion with mass in the 100 GeV range is required. It implies that the radion 
may turn out to be the lightest new particle in our model, likely to be accessible at the LHC. 
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Various studies on radion phenomenology have been performed in the past [35l [361 HOR [371 [38], I102[ 
I103[ I104[ 1105] , including recent works where radion dynamics was considered within realistic models 
of electroweak breaking with bulk SM fields [1051 11061 1107] . For example in the case of radion mass 
lighter than 2Myy, r — > 77 is a promising channel [106J, which can be also dramatically enhanced in 
the presence of Higgs-radion mixing [107] . For the case of radion mass larger than 2M\y, WW, hh, 
ZZ,ti are the dominant channels which are expected to allow for a discovery at the LHC. Thus, a 
discovery of lightish radion at the LHC and future signal at CR experiments, would yield a support 
for our class of models. 

6 Conclusions 

Indirect signals from DM annihilation in cosmic ray experiments have received a renewed attention. 
We point out that models of warped extra dimension can naturally yield a low velocity enhancement 
of the dark matter (DM) annihilation via the Sommerfeld effect. The enhancement does not rely 
on an extra dark sector, but rather is mediated via an intrinsic component of the theory, namely 
the radion with a mass at the hundred GeV range. More specifically, we studied the well motivated 
framework of a warped grand unified theory (GUT), in which the DM particle is a GUT partner of 
the top quark. Based on the Pati-Salam group, we constructed models of partial and full unification, 
which accommodate custodial symmetry protection for Z — > bb coupling. The above construction 
is consistent with electroweak precision tests for Kaluza-Klein (KK) particles with mass scale of a 
few TeV. In addition, we explore the consequences of a similar custodial symmetry protection of 
Z couplings to right-handed (RH) tau's. Such protection enables the RH tau's to be composite, 
localized near the TeV end of the extra dimension, hence having a large coupling to KK particles. 
As an aside, independently of the requirement for unification, the strong coupling between the KK 
particles and the composite tau's can lead to striking LHC signals. 

Cosmological and astrophysical aspects of our framework are discussed. We find that the dark 
matter relic abundance, as well as direct detection, constrain the viable parameter space of this class 
of models. Particularly strong constraints are found in cases where the DM particle couples to the 
neutral electroweak sector. Indirect signatures in Galactic cosmic rays (CR) are studied. We focus 
on robust observables, relatively immune to propagation model uncertainties, to test our framework. 
At present, we do not identify any clear evidence for exotic contributions. However, contrasted 
with upcoming data on the abundance of CR nuclei, near future measurements of the antiproton 
to proton flux ratio will provide a sharp probe for exotic contributions. Such contributions could 
naturally arise in our model. In case that an indirect signal is observed, measurements of the radion 
and KK particle masses at the LHC collider will provide a nontrivial test of the model. 
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A Other Pati-Salam models 

We present two other models with custodial symmetry for Zbb coupling. Just like model I (a) 
presented in the main text, both these models do not seem to fit into £0(10) representations smaller 
than 560 [62]. Moreover, even if we find a fit into a suitable larger representation of SO(W), these 
models do not have SU(5) normalization of hypercharge and hence might not maintain even the 
SM-level of unification of gauge couplings. 

A.l Model I (b): : T» R ^ and custodial for leptons 

In Tab. [HJ we first present a model with smaller SU(4) C representations than the benchmark model 
(where we had 35 of SU(4:) C ). This model has 

Y = T 3R + ^X (33) 

and hence sin 2 6' = 3/11. Also, this model has a larger value T^ R = 2, but a smaller value of 
sin#' than the model with 35 of SU (4). Such a modification tends to enhance the DM annihilation 
cross-section via Z' exchange into Zh/WW and similarly direct detection via Z exchange, whereas 
it reduces the annihilations via Z' exchange into top quarks (as per table [5]). 





SU(4) C ~ SU(3) C x 17(1)* 


SU(2) L 


SU(2) R 




10 ~ 32, 1 2 ... 
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("« t )l 


4 ~ 1 i,... 
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•3 
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2 



Table 8: Another model with custodial representations for bi and RH leptons and with non- 
vanishing v'v'Z' coupling: the subscripts denote the y/8/3 X charge. 
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A. 2 Custodial only for 6^, but not leptons 

Next, we present models with custodial representations only for bi and not for RH leptons with the 
following two motivations in mind. First of all, it is still interesting to have a scenario where RH 
leptons are not near TeV brane and thus have small couplings to Z' so that we do not need custodial 
representations for them. In this case, DM annihilates mostly into hadronic SM state, i.e., there 
is no e + signal. Moreover, we can achieve consistency with current p data, while simultaneously 
obtaining a p signal, by simply resorting to a smaller value of SE than what is used in order to 
obtain e + signal (recall that with the larger SE, the large DM annihilation into leptons was doing 
a "double-duty" of giving e + signal and maintaining consistency with present p data) . 

Moreover, there are regions of parameter space where we cannot obtain signals from DM an- 
nihilation in cosmic rays, whether e + (even with enhanced couplings of Z' to (RH) leptons) or p. 
For example, we can have heavy (~ TeV) DM as well as a heavy radion so that we do not have 
sufficient SE. Again, there is no motivation for custodial representations for RH leptons in this 
case. However, we still require custodial representations for Zbb so that it is still interesting to 
build such a unified model. 

Along these lines, the model with smallest possible representations is given in Tab. with 




and hence sin 2 9' = 3/35. 
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Table 9: Simplest model with custodial representation for bi, but not for RH charged leptons: the 
subscripts denote the -y/SS/3 X charge. 

Also, this model has T3R = 2 for v' and hence is (roughly) similar to the above model with 10 
of SU(4) as far as relic density and direct detection is concerned. 
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B The volume factor for antiproton propagation: a diffusion model 
example 



The disc+halo diffusion model for CR propagation is widely used in the literature (see e.g. 
In principle, the model allows to compute the CR densities arising from standard astrophysical 
processes on the same footing with proposed exotic contributions, such as DM annihilation. In 
practice, the model parameters are tuned on compositional CR nuclei data, which only partially 
constrains them. Here we make use of this model for two purposes: (i) to clarify some issues re- 
garding the currently fashionable "precision treatment" of exotic CR sources within a propagation- 
model dependent framework, and (ii) to illustrate the volume enhancement factor for antiprotons 
from a DM annihilation source, described in Section 14.3.41 Concerning the latter cause, we do not 
attribute particular significance to the precise numerical results, but rather consider them as order 
of magnitude estimates for the expected effect. 

We consider a cylindrical halo model with an infinitely thin disc, taking the diffusion coefficient 
as spatially constant in the propagation volume with power law energy dependence. The model 
parameters relevant in the high energy regime are L , the scale height of the cylinder, R , the 
radial extent, Dq , the normalization and 5 , the power law index of the diffusion coefficient, given 
by D(e) = D$e & . The parameters L,R,Dq,5 are constrained by B/C data, in such a way as to 
provide the measured value of the CR grammage. For relativistic energies above a few GeV/nuc, 
this constraint can be summarized as follows, 

Xesc(e) *^^g(L,R). (35) 

Above, X esc is the CR grammage, X^ sc ~ 200 pc x 1.3 m p x 1 cm -3 ~ 1.3 • 10 -3 gem -2 is the column 
density of the gaseous disc, where spallation interactions occur, and c is the speed of light. The 
dimensionless correction factor g(L,R) is given by 

it m 2R ST t f P«*\ tanh Kj) ,o R \ 
9{L,R = -rl^ Jo 2TT \ » ( 36 ) 

L fr[ v R ' wK) 

where i\ are the zeros of the Bessel function of the first kind Jo . The correction factor obeys 
g = 1 for L <C ii , and becomes smaller than one if the distance of the solar system from the radial 
edge is taken comparable to the scale height of the cylinder. For the CR grammage we adopt the 
parametrization [108 J (see also [981 1109] for earlier estimates) 

X esc = 27.5e-°- 5 gem" 2 . (37) 

We present the CR grammage in Eqs. (1351) and (1370 as a function of energy e = E/GeV . In fact, 
the grammage depends rather on magnetic rigidity, 1Z = pc/eZ . The notation is consistent as long 
as we fix our attention to relativistic antiprotons. 
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From Eqs. (f35j) and ([37|) we can deduce the following relation, 

D « 2.9 • 10- 2 (j^) e^-tyi, A) kpc 2 /Myr . (38) 

Equation (|38p may now be used in order to define sets of parameters L, R, Do, 5 , which will agree 
with high energy B/C data as long as 5 ~ 0.5 . To this end, any high energy value of e > 10 GeV 
should do. We take e = 75GeV/nuc, corresponding to the highest energy B/C measurement by 
the HEA03 mission [1101 1109] - The fact that propagation (and, in particular, diffusion) models 
must comply with the CR grammage is demonstrated, for example, by noting that Eq. (|38|) holds 
very well for the popular MIN, MED and MAX propagation models, defined in |lllj after the work 
of [115] . 

Besides the CR grammage, additional information exists on the escape time scale, found from 
measurements of radioactive CR isotopes. These data are far less accurate than the grammage 
measurements, and are given only for a limited range of energies, mostly at the ~ lOOMeV/nuc 
scale [1001 [95]. 

Different sets of values of L, R, Do, 5 , obeying Eq. ([38]) . are considered in the literature. However 
we will see that, under realistic assumptions, the diffusion coefficient does not enter into the ratio 
between the antiproton flux arising from DM and from the astrophysical background. In fact, to a 
good approximation, the only parameter which controls this ratio is the scale height of the propa- 
gation volume. We note at this point that, as the scale height L is not independently constrained, 
the DM signal to astrophysical background ratio in the disc+halo model is not constrained by the 
B/C data. We now proceed to compute the flux of antiprotons resulting from DM annihilations in 
this propagation model example. 

Neglecting losses and low energy processes and assuming steady state, the diffusion equation is 

- D(e)V 2 n = Q DM , (39) 

where n is the antiproton density. The neglect of losses kept this equation easy to analyze, at 
the price of moderate imprecision at energies below a few tens of GeV. We will return to this 
point later. Due to the homogeneity of the diffusion coefficient, the energy dependence of the 
antiproton density follows that of the source, with a trivial softening resulting from the diffusion: 
n(e, r) = e~ s f(r)Q(e, r) . We are left to deal with the spatial dependence, consisting in the function 
f(f) for which we need to derive the value in the vicinity of the solar system. 

Decomposing both n and Q in Bessel-Fourier series reduces the problem into an infinite set of 
leaky box model-like [94j equations for the coefficients. We chose a decomposition in basis functions 
which automatically satisfy the boundary conditions of vanishing CR density on the surface of the 
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cylinder. For the DM source, the decomposition reads 
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A similar decomposition holds for the antiproton density with the replacement Q m k <-> n m k . Using 



39j) we then have, for the coefficients of the antiproton density, 
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Note that the DM source is separable, 

Q-DM{r) = Qdm,p(£, r so \)n (r) , Q m k(e) = Qr>M,p(e, f so i)<? mk , 



(42) 



with (/mk the Bessel Fourier coefficients of n 2 (r). (Recall that n (r) is defined as the DM number 
density normalized to its value in the vicinity of the solar system, such that Qdm,p(c, ^soi) is just the 
local injection rate due to DM.) The antiproton density in the solar neighborhood, z = , p = r so \ , 
is thus 

aL 2 
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(44) 



Equation (|43|) allows us to obtain the specific value of the volume factor by which the DM 
annihilation source is enhanced in comparison with the production by spallation. Again neglecting 
losses, the antiproton density near the solar system, resulting from spallation, is [15] 



n 



p,spal 



PISMC 



Qspal 



P ' 



(45) 



where piSM ~ 1.3m p cm _d is the matter density on the disc. Using Eq. (|35j) and noting that 
Adisc ~ 2hpisM > where h ~ 100 pc is the half- width of the disc, we obtain the ratio of the local 
antiproton density due to DM annihilation and due to spallation throughout the Galaxy, expressed 
in terms of the local injection rates: 



%DM r QDM,p(e, ?"sol) .,, n a L 
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(46) 
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On the left panel of Figure [9] we plot the ratio of the two dimensionless correction factors a/g as 
a function of the CR halo half width L . We consider three DM halo profiles: the cored isothermal 
(ISO) and the cusped NFW, defined in Sec. 14.3,21 and the Einasto profile [113] . The ratio a/g is of 
order unity, larger for cuspy profiles compared with the cored one. To achieve faster convergence, 
we have regulated the inner cusp in the NFW and EINASTO distributions by assuming flat DM 
density for r < 200 pc. Such inner radius is not constrained by N-body simulations. We have 
verified that our results do not vary significantly as a result of increasing the regulation radius. On 
the right panel, we plot the resulting volume enhancement factor aL/gh for h = 100 pc. We find 
that, for reasonable values of L , the volume factor ranges between fy ~ 10 — 100 , depending on 
the assumed DM halo profile. 

We now comment on the neglect of losses in the discussion above. For spallation antiprotons, 
the error due to neglecting losses diminishes with increasing energy, as a result of the relatively 
rapid decrease in the grammage. For example, the error contained in Eq. (|45p due the neglect of 
losses is 25%, 10% and 5% at antiproton energies of 10, 30 and 100 GeV, respectively. Regarding 
the antiprotons from DM annihilation, the conclusion may be model dependent. However, in the 
diffusion model considered above (as well ase.g. in the leaky box model), the escape time shares 
the energy dependence of the grammage, and the conclusion is similar to the background case. In 
addition to losses by collisions with ambient matter, other low energy processes are expected to 
influence the calculation below the few tens of GeV. These phenomena include solar modulation, 
ionization losses, and even possible reacceleration or convective motion [95] . As we are dealing with 
a simplified propagation model which involves, for example, ad-hoc boundary conditions for the 
CR halo and diffusion coefficient, and an uncertain DM halo distribution, we find it useful to keep 
our expressions tractable and accurate at the high energy > 50 GeV regime, at the cost of minor 
accuracy loss below a few tens of GeV. 
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